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Abstract 

We introduce a class of discrete models for surface relaxation. By exactly 
solving the master equation which governs the microscopic dynamics of the sur- 
face, we determine the steady state of the surface and calculate its roughness. We 
will also map our model to a diffusive system of particles on a ring and reinterpret 
our results in this new setting. 
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1 Introduction 

Models far from thermodynamic equilibrium, may not obey detailed balance and 
hence their steady state can not be described by Gibbs Measure, and one often needs 
to study closely the stochastic dynamics and obtain the steady state directly by 
simulation, numerical or exact solution of the microscopic dynamical equations. (For 
recent reviews see |2| and references therein). The microscopic dynamics which is 
usually encoded in a set of transition rates between different configurations, usually 
bring about strong long-range correlations in one dimensional systems. This is the 
case, for example in driven diffusive systems the prototype of which is the Asymmetric 
Simple Exclusion Process|Q. In such cases then, the mean field solution is not exact 
or else is exact only under certain conditions on the rates. 
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In recent years, much attention has been paid to the exact solution of stochastic 
processes on one dimensional lattices. Among the methods developed, the Matrix 
Product Ansatz (MPA), has proved particularly fruitful, since it has both led to 
simple and nice solutions of the previously known models Q and has also led to 
various generalizations [5-12]. 

The Matrix Product Ansatz is best understood within the operator formalism for 
Markov processes on one dimensional chains in which the steady state of the process, 
turns out to be the ground state of a suitable Hamiltonian with nearest neighbor 
interactions 

L 

H = Y^ hk,k+i. (1) 

k=l 

According to MPA, the steady state probability P(ti,T2, ■ ■ ■ tj^) for the configuration 
{ti,T2, • • • tl) is written as the trace of a product of noncommuting matrices: 

Piri,T2, ■■■tl) = ^^tr (^D{n)D{T2) • • • D{tl)^ (2) 

(Here we are taking a closed system with periodic boundary conditions.) It is then 
easy to show that the above probabilities are the steady state ones, provided that the 



matrices D{t) satisfy a set of relations encoded compactly as:|13| 



/iDoD = X®D-D®X (3) 

Here D and X are matrix valued column vectors. 

The solution of these equations are however far from trivial and so far relatively 
few systems have lend themselves to this type of solution. In cases where we take 
D(t) to be one dimensional and accordingly denote them by d(r), we are in fact 
searching for exact mean field solutions. Even in this case, (^) is a set of coupled 
nonlinear algebraic equations and are generally difficult to solve, or else so many con- 
straints of consistency are imposed on the rates to render the process uninteresting 
and non-physical. 

One can however reverse the problem [14, ^, 1^, 15] and search among quadratic alge- 



bras for those classes that may correspond through the MPA formalism to a stochastic 
Hamiltonian. 

Mathematically, a through search and classification is quite difficult, and one is better 
to start from a restricted class of algebras or processes which seem to lead to physi- 
cally interesting models. 

( For a partial classification of algebras, which pertain to diffusion and exchange of 



particles see [|15|). 

In this paper we take all the matrices to be one dimensional which means we are 
considering those processes for which mean field solutions are exact. As far as the 
number of particles of each species on a ring is conserved, a one dimensional repre- 
sentation leads to a rather uninteresting steady state, namely, one in which all the 
configurations have equal weights. However if there is no such conservation law, then 
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these mean field solution may be quite nontrivial and contain interesting physics, as 
we shall see. 

Our starting point is to write a set of microscopic dynamical rules pertaining to 
a generalization of the so-called zero-range contact process |3[(see also [16] for more 
references) and ask if the steady state of the corresponding stochastic process allows 
an exact mean field solution (i.e: a one dimensional MPA solution). We find that an 
interesting and rather general class of these processes allow such solutions, namely 
those processes in which the rates of transitions between neighboring sites can be 
factorized in two pieces each pertaining to one of the sites . These processes can 
have two interpretations, one in terms of surface relaxation and the other in terms 
of diffusive systems of particles. For definiteness we adhere to the first interpretation 
and determine exactly the steady state of the surface and calculate the roughness of 
the surface. We study in detail two members of the class. 



2 Models for surface relaxation 

We consider a ring of L sites, each site of which can accommodate a number of 
particles on top of each other. At most p particles can be accommodated in each site. 
To each site k of the lattice {k = 1, • • • L), we assign a stochastic variable which 
shows the number of particles in that site, i.e. (r^ = 0, 1, • • - p). 

The microscopic stochastic dynamics which governs the relaxation of the surface 
are described according to the following transition rules between the stochastic vari- 
ables pertaining to any two adjacent sites k, and k + l(Fig. 1): 



i j — > i — 1 j + 1 with rate fiij, (4) 
j i — > j + 1 i — l with rate fiij, (5) 

where i and j represent the values of r^,, and t/^^i. Note that the model has left-right 
symmetry and we are taking time to be continuous, hence we speak of rates rather 
than probabilities. Obviously we should have the following constraints on the rates, 
namely fxoj = fijp = Vj. Moreover we require that the higher the difference 
between the heights of adjacent columns, the more probable is the above transition. 
Thus we require that fiij be an increasing function of {i — j). 

Note that the surface does not grow on the average and the total number of particles 
is constant which we take to be equal to N. If the relaxation of the surface is much 
faster than the rate of fall of particles , then this model can also describe the growth 
of the surface, in this case time will be proportional to the number of particles N. 



3 Interpretation as a diffusive system of particles 

Our model can also be interpreted as a model for particles hopping on a ring. One 
can now think of L particles on a ring of L + sites, where each particle k has a 
number of vacant sites ahead. In this view, the transition rules are interpreted 
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as follows: A particle whose distances to its left and (right) neighboring particles are 
i and (j), hops forward with rate and backward with rate ^ij 

The increasing property of /i, j in terms of (i — j) is still plausible in this setting, 
specially if we think of traffic flow, since a car in the above condition accelerates if 
i — j is small and decelerates if i — j is large. In this interpretation the probabilites (2) 
give the distribution of headways and the roughness designates how uniform the flow 
is. For reviews on traffic flow and (particularly its statistical mechanical aspects), see 
[|l7|], and (|]l^]) respectively. 



The above process can now be studied as an interacting stochastic system. For 
this study we use the operator formalism of Markov processes, where to each configu- 
ration {ti,T2, ■ • • Tl) of the system a state vector \ti,T2, ■ • • tl > is assigned and each 
probability distribution P, is encoded in a vector 

\P>:= P(.ri,T2,---rL)\Tur2,---TL)> . (6) 

The totality of possible state vectors span a complex Hilbert space, denoted hy Tl := 
^(jp+i^®L^ where C^+^ is the complex space spanned by = { |0), |1) ... \p)}. The 
stochastic Hamiltonian describing the above dynamical rules, is given by 

L 

W = X! ^k,k+i, (J) 

k=l 

where 

~ \i-lj + j\ + \j + li- i\ - \i j){i j\ - \j i){j i\ 1 . (8) 



Inserting (P) with the choice D = column{do,di, ■ ■ ■ dp) and X = into equation 
(^) we obtain : 

fJ'i+i,j-idi+idj-i + fij+i,i-idi-idj+i = {fiij + fij^i)didj (9) 
A sufficient condition for solving these equation is the following: 

fJ-j+i^.-idi-idj+i = Hijdidj \fi,j, (10) 

For general rates these relations may not have any non-zero solutions, however if 
we take the rates to be of factorized form, namely 

fiij = HiXj, (11) 

where iJ,i and Xj are arbitrary non-negative real numbers, then the set of equations (pX|) 
have the following solution for dj's: 

J 1 J -^0 , ^1-^0 J A2A1A0 Ap_i---AiAo , , 
do = 1, di = — , d2 = 03 = ••• dp = (12) 
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Remark: 

The symmetry of the process is crucial in the above derivation, for a totally asym- 
metric process, the rates will be highly constrained. 



Note that the factorization property is the only requirement that exact solvability 
imposes on this model. We are free to choose the dependence of fii and Aj on i as we 
wish. It is for physical reasons that we require that 

/Up > /Up-i • • • /^i > A^o = (13) 
Ao > Ai > • • • Ap-i > Ap = (14) 

In this way the rates /ijj will be increasing functions of (i — j). Furthermore the 
probability of a particle jumping from a lower height to a neighboring higher position 
will be quite low. We can now calculate the probability of all configurations. From 
(0) we find: 

Piri,T2, ■■■Tl) = J^^^din)d{T2) ■ ■ ■ d{TL) (15) 

Where Z{L,N), the partition function, is given by: 

Z{L,N):= d{Ti)d{T2) ■ ■ ■ d{TL) (16) 

Due to the commutativity of (i(r)'s, the probabilities depend only on number of 
columns Ki of a given height i, and not on the position of columns. Thus ([T5| ) can 
be rewritten as: 

^(^0. K,,-.K,) = . . . (17) 

where the combinatorial term account for various possible places of the columns, and 
(16) can be re-expressed as: 

Where by " we mean a sum subject to X^iLo ~ ^^=o ~ ^- Our 

analysis will be simplified if we consider a grand partition function Zg{L, s): 

pL 

Zg{L,s) = 5]5^^(L,iV) (19) 

N=0 

= ido + dis + ...dps'ff =: iz{s))^. (21) 
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The parameter s plays the role of fugacity. The average density of particles n := ^ 

and the average ratio of columns of height i, that is ; ki := are given respectively 
by: 



s d ^ _ . d ^ , . 
n = -— In Zg{L,s) = s— Inz(s) 



h = ^^lnZG{L,s) = ^ (22) 
L odi z[s) 



The roughness of a surface is usually defined as |19|: 



n 



(23) 



where h(x) is the height of the surface at point x. This can be rewritten as: 

w 



\ 



1 



L 



L 



A simple calculation shows that: 

72 



ds ds 



(24) 



(25) 



We now consider two examples which lend themselves to elementary analysis. 
3.1 Example A 

In this first example, we take the rates to be as follows: 

Hi = i Xi= p-i (26) 

The above assignments satisfy all the requirements put forward on physical grounds. 
For example the rate of fall of a particle from height p to height is p^, while the 
rate of jump of a particle from a height 1 to p — 1 is 1. Moreover /io = Ap = 0. In 



this case we find from (12): 



di 



p\ 



i\{p — i)! 



and 



z s 



1^^' r^S 



From (|2^) we find the average number of particles: 

ps 



(1 + sY 



n 



n 



or 



1 + s p — n 

Also from (|2^) we find the fraction of sites with height i to be: 



ki 



P 

i 



[l + s)P 



(27) 
(28) 

(29) 
(30) 
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Rewriting this in terms of the number of particles n, we find: 

k^ = ^(^''^y{p-nr-\ (31) 

which implies a binomial distribution for the columns of different heights. 
The roughness is calculated from (p5|) and (|2^) to be: 

W = ^n(l-^), (32) 

If the relaxation is much faster than the rate of fall of particles, we can take n ~ t and 
p ~ T where T is the time for growing a layer of uniform thickness p, then the above 
equation implies, that at the early stages of the growth when n <^ p, the roughness 
increases as ~ tz, in the middle stages of the growth, n is comparable with p, it 
increases linearly with time, and finally at the final stages when the number of holes 
p — n becomes a good measure of time T — t, the roughness again decreases as the 
square root of T — t. 



3.2 Example B: 

In this example we consider the following rates: 

fii = 1 yi ^ 0, Xi = 1, yi p fiQ = Xp = (33) 

That is, the rates do not depend on the height of source and target columns. We find 
from (^) that do = di = ■ ■ ■ dp = 1 and thus: 

1 _ oP+i 

z(.) = ^— — (34) 

from (^) we find: 

Large p-limit Some insight is obtained if we look at the large p limit of the roughness 
for the two cases s < 1 and s > 1. 

When s < 1, we find from (|3^), that z(s) ~ and using ( p2| ) and ( p5| ) we obtain: 



n~ Wc^^ — =Jn(l + n) (36) 

I — s (1 — s) * 

On the other hand when s > 1, we find that Inz(s) ~ (p+ 1) In s — ln(s — 1), which 
when inserted into (22) and (pS]), yields respectively 

n~p and 1^ ~ = J (p — n)(p — n + 1) (37) 

s — 1 s — 1 * 
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Moreover from (35) we find that in the large p hmit, when s < 1 or n < p/2, 



and when s > 1 or n > p/2 



71 

(l + n)^+i ^^^^ 



' (39) 



(l+p - n) 

It is interesting to note that despite the very different nature of the transition rates in 
models A and B, in both models, the formulas obtained for the roughness imply that 
the roughness first increases as t^/^ and then increases linearly and finally decreases 
at the final stage again as (T — t)^!"^ . 

These two regimes are separated by the case when s = 1, or n = |, where we find 
for all arbitrary p's, ki = -pp^^ meaning that all heights are equally probable. In this 
case the roughness turns out from (25) to be: 




(40) 

Fig. (3) shows the value of roughness in terms of the fugacity for p = 10. It is seen 
that there is a sharp peak in the roughness at n = |, when half of the capacity of the 
lattice is filled with particles. 

All these hint that the point s = lorn = |isa special point, possibly a point of 
phase transition. 

Acknowledgement: V. K. would hke to thank V. Rittenberg for his valuable 
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the manuscript. 
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Figure 1: In the surface relaxation model, a particle drops with rate 



10 




Fi gure 2: In the diffusive system, a particle hops with rate Hij 
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Figure 3: The square of roughness for p=10 
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